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Abstract 

We explore the possibility that a new physics interaction can pro- 
vide an explanation for the knee just above 10 6 GeV in the cosmic 
ray spectrum. We model the new physics modifications to the total 
proton-proton cross section with an incoherent term that allows for 
missing energy above the scale of new physics. We add the constraint 
that the new physics must also be consistent with published pp cross 
section measurements, using cosmic ray observations, an order of mag- 
nitude and more above the knee. We find that the rise in cross section 
required at energies above the knee is radical. The increase in cross 
section suggests that it may be more appropriate to treat the scatter- 
ing process in the black disc limit at such high energies. In this case 
there may be no clean separation between the standard model and new 
physics contributions to the total cross section. We model the missing 
energy in this limit and find a good fit to the Tibet III cosmic ray flux 
data. We comment on testing the new physics proposal for the cosmic 
ray knee at the Large Hadron Collider. 



1 Introduction 

A power law energy dependence going approximately as E~ 2 ^ characterizes 
the high energy cosmic ray spectrum from 10 GeV to 10 6 GeV. In the range 
10 6 GeV to 4 x 10 6 GeV, the energy dependence changes to E~ 3 , continuing 
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with this slope to 10 GeV or so. The downward bend just above 10 GeV 
is called the "knee" in the spectrum. This knee phenomenon in the cosmic 
ray spectrum [I], observed by many experiments over many years [21 13 HI [5] , 
and reviewed recently in jS], still lacks a convincing explanation. It is gener- 
ally believed to be of astrophysical origin, caused by leakage of cosmic rays 
from the galaxy or due to an upper limit on the energy that cosmic ray 
particles can attain at supernova remnants, the galactic acceleration sites 
[BJ. The center of mass energy corresponding to the knee is several TeV 
in the proton-proton (pp) system, just above the highest energy laboratory 
measurements of pp cross section and 2 orders of magnitude above the high- 
est pp measurements, but close to where one may argue that new physics 
such as supersymmetry, technicolor, or low scale gravity may begin to make 
contributions to pp scattering. Without adopting any particular symmetry 
or dynamics, one may thus speculate that some change in the fundamen- 
tal hadron interactions might be the cause of the knee in the cosmic ray 
spectrum[71 13 QUI QH E2]. To explain the existence of the knee, one may 
propose that the incident spectrum obeys a single power law, but beyond a 
certain energy the collision between the cosmic ray primary and the target 
atmospheric nucleus is dominated by interactions that lead to an enhanced 
cross section and production of unobservable, weakly interacting particles, 
leading to missing energy. Hence the energy deduced from the shower does 
not correctly estimate the energy of the primary, leading to the knee in the 
spectrum. 

In the present paper we pursue the consequences of this hypothesis. In 
particular we are interested in determining whether the measured nucleon- 
nucleon cross section at ultra high energy using cosmic ray data |13[ [14] is 
consistent with the cross section required by this mechanism. Furthermore 
we are interested in determining how the new physics explanation of the 
cosmic ray knee may be tested by observations at the Large Hadron Collider 
|15j . Analysis of the impact of changes required in the nucleon-nucleon 
cross section within this framework, analysis necessary for comparison both 
cosmic ray and collider data, is so far not available in the literature. We 
present here a step toward filling this gap. 

We first assume that the cross section for new particle production adds 
incoherently to the standard QCD nucleon-nucleon cross section, with a 
large portion that does not contribute to the observable shower. The new 
particle production process turns on only at very high energy and hence 
this would predict a sudden increase in the effective nucleon-nucleon cross 
section beyond the threshold for new physics. This simple model for the cross 
section, however, may not capture all the essential physics of the process. At 



2 



high energy the nucleon-nucleon cross section becomes sufficiently large so 
that it may be more appropriate to apply a black disc model to compute the 
proton-air nucleus cross section, relevant for cosmic ray collisions. In this 
case the proton-air nucleus total cross section has very weak dependence 
on the pp cross section. The scattering involves multiple interactions of 
the incident proton with the nucleons in the nucleus, which itself involve 
multiple interactions of partons. If there exists a new physics interaction 
with an interaction strength comparable to the QCD interaction at energy 
of the order of a GeV, it is very likely that almost every proton-nucleus 
collision involves at least one new physics interaction. This will imply a 
considerable modification of the first approach outlined above. 

An additional complication is that the nucleon-nucleon cross section is 
dominated by soft hadronic physics and hence cannot be computed from first 
principles. Although a precise first principle calculation of the cross section 
may not really be required for a reliable analysis of the cosmic ray data, one 
does require a good understanding of the particle production rate in order 
to accurately extract the energy of the primary particle from observations. 
In the absence of a first principles calculation, one must extrapolate a model 
calculation to high energies. Such an extrapolation may lead to large sys- 
tematic errors in the extraction of basic parameters such as the energy of 
the primary particles. 

We adopt several simple models for the new physics cross section and 
explore a wide range of parameter values in each. These cross section models 
are given in Section 2. We use the recently published, high statistics data 
for the cosmic ray spectrum from the Tibet III experiment [5] to implement 
our program. Our fits to the cosmic ray spectrum in the knee region and 
the resulting prediction of the pp cross sections are discussed in Section 

3. The new physics scale and the fraction of missing energy must be such 
that the Tibet data are well fit and at the same time, reasonable agreement 
with the published ultra high energy cross section data [131 E] is obtained. 
In every case, we find parameter sets that give a good description of the 
Tibet III spectrum with a single power law injection spectrum. However 
the required cross section values are found to be much larger in comparison 
to the values extracted from cosmic ray data. This is discussed in Section 

4. In Section 5 we discuss the black disc limit of the cross section. Here the 
multi-particle scattering obscures the link between the p-air nucleus total 
cross section and the pp total cross section. In this limit it is not possible to 
write the total cross section as the sum of the contribution due to standard 
model and that due to new physics. We make a fit to the cosmic ray flux by 
extrapolating our fit to the accelerator pp total cross section to the energy 
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range of the Tibet III data, modeling the fraction of energy loss with a power 
law dependence on energy. We find a good description of the Tibet III data 
with this simple, alternative picture as well. Finally we conclude in Section 



2 The total pp cross section: old and new physics 

We assume that the total proton-proton cross section derives contribution 
from both known physics (Standard Model, laboratory data) as well as some 
new physics effects. We express the total pp-cross section, a tot as, 



where a sm is the pp-cross section calculated using known physics while a np 
is the cross section due to new physics. 

This model is analogous to those used in early studies of the growth of the 
pp and pp cross sections with energy [161 113 (HI 021 1201 EI] ■ There the cross 
section was represented as atot = &soft + a QCD, where <J S oft incorporates 
the long distance, non-perturbative component and ctqcd incorporates the 
increasing number of perturbatively calculable jets. This early picture treats 
the "soft" part as a constant, fixed at energies just above the resonance 
region, and ascribes the growth to the sum of perturbative QCD "mini- 
jets". The picture, though simple and intuitively appealing, basically the 
one adopted in [TJ, runs into contradiction with partial wave unitarity [22], 
which can be addressed by adopting a diffractive approach, implemented 
by an impact-parameter representation [23U24"]; later studies incorporated 
unitarity constraints [25l ES E3 Ell E3 IHofT We will explore two models 
here, where the "soft" term will be replaced by several parameterizations of 
the "standard model" part of the high energy cross section and the "mini- 
jet" portion of the cross section by models that produce a rapid rise of the 
cross section in the PeV region. 

If the new physics process involves production of new, exotic particles 
that interact weakly with the atmospheric environment, they will not con- 
tribute to shower energy. For incoming cosmic rays with energies above the 
new physics scale, the observed shower energies will then be lower than the 
energies of the cosmic rays that initiated them. There would be a depletion 
of events above this scale, and a corresponding enhancement below. There 
would be an apparent "knee" in the cosmic ray spectrum. 

1 Refs. 22 - 26 incorporate pp and pp data as low as 6 GeV and also discuss simultaneous 
fits to the cosmic ray cross section data. 
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In Eq. Q] we display our working model for atot ■ This model assumes 
that one can cleanly separate the events in which either only the standard 
model interactions contribute from those where only the new interactions 
contribute. This leads to a convenient probabilistic approach that we outline 
next. 

At high energies a cosmic ray proton primary with energy E interacts 
with the atmospheric nuclei with probability 

(E)/(a sm (E) + a np (E)) (2) 
through the standard channel and with probability 

P np (E) = a np (E)/(a sm (E) + a np {E)) (3) 

through the new channel. Here the a sm and a np refer to the cross sections of 
primary cosmic ray particles with air nuclei. However we shall assume that 
we can replace them simply by the corresponding values for the pp cross 
section. This is justified as long as o sm and o np scale proportionately as we 
go from pp to p-nucleus or nucleus-nucleus scattering for light nuclei. The p- 
nucleus cross section can be computed using the Glauber multiple diffraction 
theory J3TJ [32], |27j . One finds that the p-air cross section, computed in Ref. 
|27j . depends almost linearly on the pp cross section as long as the cross 
section is not very large, suggesting that it may be reasonable to apply the 
Glauber theory to the two cross section a sm and o np separately and to adopt 
the scaling proportionality. However, as discussed in the introduction, this 
assumption may break down at ultra high energies, where the cross sections 
become very large. 

We next assume that if the cosmic ray primary interacts through the 
new channel then total energy detected is E = yE', where E' is the total 
energy of the incident cosmic ray and y < 1 is an energy loss factor that 
characterizes the new physics, assumed here to vary slowly enough to be 
effectively energy independent. Because proton collision can be initiated by 
either the new physics or the standard model interaction, we find it useful 
for later application to introduce an effective, energy dependent y e // 

y eff (E) = P sm (E) + yP np (E), (4) 

where y is the value that follows from the fit of the model data, and y = 1 
for the Standard Model. Clearly for E < E np , where E np is the threshold 
scale of new physics, y e fj = 1. 

Next let 4>{E) = ME -1 exp(—E/E cut ) be the incident galactic cosmic ray 
flux, with 7 « 2.7 and the cutoff E cut ~ 10 8 to 10 9 GeV, representing the end 
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of galactic sources to the CR spectrum. Here M is an overall normalization. 
The observed intensity at energy E is then given by the incident intensity 
at E times the probability that the observed showers are initiated by a 
standard physics collision at that energy plus the incident intensity at energy 
E/y times the probability that the observed showers are initiated by the 
new physics interaction at the energy E/y, which takes into account the 
energy loss due to new physics. We can express the observed differential flux 
spectrum, (p b s (E), as 

/•oo 

cf> obs (E) = <t>{E)P sm + / dE'6(E - yE')0(E>)P np (E>). (5) 

Pulling out an overall factor of the incident intensity at E, we write the 
observed flux at E as 



tobs(E) = NET* eM-E/E cut ) [P sm (E) + y^ e -(E/yE cut )(l-y) p np{E/y ) . 

(6) 

With this general framework, we are equiped to determine the new physics 
parameters required to fit the data for the cosmic ray spectrum covering the 
knee region. To proceed, we need to specify the laboratory-measured pp and 
pp total cross sections and the generic, new-physics models we will use for 
our study. 



2.1 The laboratory-measured, "standard model" pp cross 
sections 

Let us first fit the proton-proton cross section data above 50 GeV [33J from 
accelerator measurements, but excluding the high energy points reported by 
the Fly's Eye [13] and AGASA [14], which lie in the energy range 2 x 10 7 
GeV < E < 4.8 x 10 8 GeV, an order of magnitude and more above the 
energy where the knee in the spectrum appears at 2 — 4 x 10 6 GeV. We can 
then ask whether the extrapolation of the fit to the laboratory cross section 
measurements plus the new physics cross section found from adjusting its 
parameters to fit the CR spectrum with an assumed flux that falls with a 
fixed power law, namely E~ 2m below the knee, according to the Tibet III 
data, reproduces the direct CR cross section measurements lying above the 
knee region of the spectrum. 

In order to have a reasonable fit to the laboratory energy data [33], we 
assume the following "Froissart bound saturation" form for o~ sm (s) [281 [291 
[30] . known to yield good fits to the laboratory data above about 10 GeV. 
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Data set Cq C\ C2 X 2 


d.o.f. 


pp > 50 GeV 46.60 -6.875 1.472 10.5 


29 


pp > 50 GeV +pp 41.19 -3.863 1.077 15.6 


32 



Table 1: First row: parameters, \ 2 an d number of degrees of freedom for 
Log 2 {E I Eq) fit, Eq. 7, to fixed target pp total cross section data with E > 
50 GeV. Second row: same data plus 3 collider points. 

The fit function is 

osm = C + Ci log(s) + C 2 log 2 (s) , (7) 

where s = 2m 2 + 2m p E is the center of momentum frame total proton- 
proton energy squared. Since we are strictly interested in the high energy 
regime and a fit to the cosmic ray cross section data we fit only the labo- 
ratory data set above 50 GeV. The resulting parameter values are given in 
Table [TJ As an additional exercise, we include three pp points from collider 
measurements, which lie several orders of magnitude above the available pp 
laboratory values. Including these points that extend nearly to the energy 
corresponding to that of the knee, we find the parameter values given in the 
second row of Table [TJ 

The plots of the cross section fits versus data are shown in Fig. [TJ which 
includes the highest energy points from the cosmic ray data, not included 
in the fit. 

2.2 Two new physics cross section models 

We start our discussion with a generic new physics cross section, which 
asymptotes to a constant times the extrapolated standard model cross sec- 
tion cr s?Tt We chose a simple parametrization that may be expressed 
as. 

a np (E) = aa sm (E)T(E)g(T{E)) (8) 

Here a is a dimensionless constant, t(E) = Mq/s, s ~ 2m p E and Mq 
is the scale of new physics. In the nucleon-nucleon collision, the factor 
t{E)q{t(E)) represents the differential gluon-gluon luminosity having in- 
variant mass-squared s > Mq |34j . The function g(r(E)) takes the form 

g(r)= ! dxf(x)f(r/x), (9) 
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Figure 1: pp cross section above 50 GeV plotted against energy showing fits 
(see text) that include three collider pp data points, lower curve, and do not 
include them, upper curve. The cosmic ray data are displayed, but are not 
included in the fit. 



The function f(x) is parametrized as |35j 

f(x) = v^+T^ — (10) 

x 

In the limit E/(Mq /2m p ) — > oo, r — > but Tg[r) — > 1, so rg(r) acts as a 
"step" function that rises from at energies below the new physics threshold 
at E = Mq/(2top) to 1 at higher energies. In this limit, a np (E)/a srn (E) — >• 
a; for this model the total cross section approaches (1 + a)a sm at energies 
much larger than Mq/(2771 p )H 

Our second model is a generic parameterization that copies the \og{E) 
plus log 2 (£) fit that is very successful at describing the laboratory total pp 
cross section, as we see explicitly in the following section. Introducing a new 
scale E n , we assume the new physics total cross section to be of the form 

a np (E) = Dt \og(E/E n ) + D 2 log 2 (E/E n ), (11) 

with a np (E) = if E < E n , and we fit the flux to the form of Eq. 6. The fit 
parameters are E n , the new physics scale, y, the energy loss parameter, A/", 

2 This is the same log 2 (E) asymptotic behavior as the "Froissart bound" fit to the data 
[25] . but with a strength 1+a times as great. 
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the normalization, D% and the coefficients of the linear and quadratic 
\og{E / E n ) terms in the new physics cross section. 



3 New physics models fit to Tibet III data 

First we explore the parameters of the model by fitting it to the average of 
the Tibet III QGSJET-HD [361 EZ] and SIBYLL [38] data in the knee region, 
using the difference of each from the average to estimate a systematic error. 
The total errors are then taken to be the rms value of this systematic error 
estimate and the quoted, essentially common, statistical errors associated 
with the two analysis packages. 

We compare model 1 with the data by fitting the form 



Vobs 



(E) = NE-* exp(-E/E cut ) 



a yl-l e -E/yE cu t+E/E cut 



1 + aTg(r) a + l/[(ry)g(ry)] 

(12) 

to the data, finding values for parameters J\f, a, Mo, y, and E cut . Because 
of our definition of model 1, cr sm does not appear in (f) bs- An example fit 
and comparison of the model with the data is shown in the dashed curve in 
Fig. [2j Here, following [7] for illustration, we set the parameters n = 6 and 
we set 7 = 2.666, the value we obtain for the region E < 10 6 GeV. 

The fit is good, creating the features of the data with a uniformly falling 
power law spectrum cut off at 10 10 GeV, with a new physics cross section 
that is proportional to the SM cross section with an energy dependent coeffi- 
cient added to the SM cross section. The fit finds a new physics characteristic 
scale just below 900 GeV and the energy loss parameter of y ~ 0.2, which 
ensures that most of the energy from a cosmic ray collision is invisible. 

Next we also fit the data using 



IJ obs 



(E) = A/\E~ 7 exp(—E/ E cut ) 



R{E/y)y^ l e- E /y E ^ +E / E - 



1 + R{E) 1 + R{E/y) 

(13)" 

where R(E) = [D x \og{E/E n ) + D 2 log 2 (E/E n )]/[C + d log(s) + C 2 log 2 (s)] 
and s = 2m 2 + 2m n E, that follows from model 2. In Fig. [21 the solid curve 
shows the case where only the pp data are used for the laboratory physics 
input, which are fit separately with the form a^ 1 (E) = Cq + C\ log(s) + 
C2 log 2 (s), as described in section 2.1. When we include collider data points 
for pp [33] in the SM fit, the result differs little from the result for only pp 
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Figure 2: Energy spectrum with E along the x-axis and E 2 ' 5 x flux along 
y. The dashed curve is the model 1 fit and the solid curve is the model 2 fit, 
using only pp laboratory data as input, with no constraints from the cosmic 
ray cross section estimates. 

data shown in Fig. [2l For the fit without the collider points, the essential 
parameter values are E n = 2.17 x 10 6 GeV and y= 0.20, while for the fit 
with the these points, we find E n = 3.04 x 10 6 GeV and y = 0.20. Again 
the fits are quite good. 

4 Fitting the all proton-proton cross section data 
and Tibet III data with new physics 

Referring to Fig. [TJ we see a gap between the highest energy pp lab data 
and the lowest energy CR data, a gap which contains the energy range of 
the knee in the CR spectrum. The question we address next is: "What is 
the behavior of the total pp cross section, new physics plus measured lab 
values, implied by our new physics analysis of the Tibet III spectrum, which 
spans the energy range from E = 1.12 x 10 5 to E = 1.78 x 10 8 ?" This range 
includes the collider cross section measurements and most of the CR cross 
section measurements. 

In order to compare the cross sections extracted from cosmic ray data 
to our predicted values, including new physics contributions, we must ap- 
propriately rescale the energy assigned to these cross sections values. This 
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is because the analyses [131 E] that led to the highest energy cross section 
points in Fig. [3] assign primary, collision energies to the events based on 
standard analysis of totally inelastic p-air collisions that produce showers 
containing all of the incoming primary cosmic ray energy. If instead, there 
is substantial energy missing, energy in a form not accessible to the detec- 
tors, then the primary energy must be higher than contained in the observed 
shower. This implies that the energy assigned to events is too low by the 
factor 1 / y e f f , where y e f / is the effective energy loss parameter as defined in 
Eq. HI Therefore the energy of each cosmic ray cross section point must be 
rescaled by this factor. 

In Fig. O we show the total cross sections, new physics plus extrapo- 
lated laboratory fits, plotted versus energy for new physics model 2. This 
is the case where the fitting function is linear plus quadratic in Log(E/.E n ), 
as described in the preceding section. We show both of the cases we con- 
sidered, the pp lab data only and the pp lab data plus three pp points. The 
details differ, but in both cases the cross section required to produce the 
apparent knee in the CR spectrum rises explosively at energies below those 
where cosmic ray measurements constrain the cross section, albeit with large 
uncertainties. The model 1 case is quite similar to that shown for model 2 
in Fig. El 

We find that the predicted cross sections show dramatic disagreement 
with those measured from cosmic ray data. Hence the simple cross section 
model, Eq. (H is ruled out by the cosmic ray cross section data. However 
as discussed in the introduction, when the cross section values become so 
large, we expect that a typical pp collision involves multiple scatterings of the 
elementary quarks and gluons. Hence once the energy becomes sufficiently 
high a typical process might involve many standard model interactions at 
the parton level and perhaps one (or few) new physics interactions. In this 
case there may be no clean separation between the standard model and new 
physics contributions to the pp cross section. This distinction gets further 
clouded when we consider the proton-nucleus collisions, relevant for cosmic 
ray showers. We take up these questions in the next section. 

5 Black Disc Limit 

At ultra high energies the simple formula, Eq. [H employed so far for pp scat- 
tering cross section may no longer be applicable. As discussed in Section 
1, in this limit it may not be possible to cleanly separate the contributions 
entirely due to standard model interactions from those due to new physics 
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Figure 3: pp cross section data above 50 GeV compared to model 2 cross 
sections fitted to laboratory data with 3 collider pp data points, lower curve, 
and without these points, upper curve. The new physics parameters are 
determined by fitting model 2 to the Tibet III flux data. The high energy 
data points show the rescaled (see text) cross sections extracted from the 
cosmic ray data. The red and blue points employ rescaling corresponding 
the upper and lower curve respectively. 

interactions. Lets assume that the new physics interaction strength is com- 
parable to that of the QCD coupling at the energy scale of a few GeV. In 
this case a typical pp collision may involve many standard model interac- 
tions along with some new physics interactions. This is even more true in 
the case of proton-nucleus collisions relevant for cosmic rays. In this case, as 
the cross section grows the incident proton undergoes multiple scatterings 
with the target nucleons such that the proton-nucleus cross section becomes 
nearly independent of the underlying pp cross section above some black disc 
threshold |39|, 132] . Thus once we cross the threshold for production of the 
"invisible" particles associated with new physics it may be more natural 
to assume that almost every collision involves a certain amount of energy 
loss due to new physics interaction. The amount of energy loss per collision 
would increase with energy due to increase in probability for new physics to 
contribute. 

To implement the idea just described, we assume that the process is 
described by a single total cross section cr(E), which has a threshold for 
production of new invisible particles and, in general, a change in magnitude. 
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The energy loss parameter, y, is taken to be a non-trivial function of energy, 
f(E), rather than a constant as in the earlier discussion. We will take it 
to be a power of energy shortly, but consider it to be unspecified for the 
moment. The number of events per unit of energy, time, solid angle and 
area, suppressing acceptance, can be written 



4>obs{E) 



dE' 



6{E T - E')5{E - E')(j)(E') 



+6{E' -E T )5(E- f{E'))<j>{E') 



(14) 



where Et is the threshold for onset of new physics effects. Performing the 
trivial integration, we find 



4(E) obs = 4(E)0(E T -E)+ fl E °™ d(E (E) - E T ), (15) 

df{E (E))/dE 

where Eq is the solution to E = f(E'). Assuming a power law form f(E) = 
E 1+P /E T , we find a solution for the power needed to produce the break at 
the knee in the Tibet III spectrum with p = —0.24, and Et = 3.68 x 10 6 
GeV. Within the uncertainties in these extrapolations i 101 , there is room for 
a new physics contribution to the pp cross section of the order of a factor 
two or so. If it has an energy -dependent missing energy fraction that grows 
with energy as modeled here, it can plausibly produce the observed knee in 
the cosmic ray spectrum, as illustrated for the case of a simple power law 
dependence, in FigHl 

These results show that the situation changes when the cross section 
is saturating the multiple scattering limit, where the cross section does not 
cleanly separate into several pieces at each energy, and the new physics frac- 
tion of visible energy increases with energy. In this case we are unable to 
predict the dependence of cross section on energy. For example, we could 
change the cross section, while at the same time modifying the form of the 
energy dependence of the missing energy. Hence we find that we are unable 
to test the new physics hypothesis to explain the cosmic ray knee simply by 
total cross section measurement at the LHC. A more detailed analysis which 
involves studying the interaction strengths of new physics processes is re- 
quired. First, a systematic failure to match the total final state energy with 
the LHC collision energy would have to be established to support the pic- 
ture we present here. Then, taking the TOTEM experiment for definiteness 
1151 . the measurements of the numbers of elastic and inelastic events, N e i 
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Figure 4: The Tibet III data is fitted with a cross section model that is 
simply an extrapolation of the lower energy fit up through the pp collider 
data. At a threshold Et, new physics sets in that produces a fraction of 
missing energy that is controlled by a power p. The threshold is at 3.68 x 10 6 , 
and the p value is —0.24, in the curve through the data shown in the figure. 

and Ni ne i, along with dN e i/d\t\\\ t \ = o, would have to be corrected to account 
for the missing energy to find the total cross section via the optical theorem 
and to determine the luminosity. Conversely, if the cross section and lumi- 
nosity are determined by these measured quantities without correction, then 
if an independent luminosity measurement is inconsistent with that deter- 
mined by N e i and N ine i it could be evidence for new interactions producing 
undetected energetic particles, i.e. missing energy. Such an independent 
luminosity determination uses the precise measurement of the horizontal 
vertex distribution at the interaction point (IP5 of LHC), as mentioned in 
the conclusions of the second paper cited in |15| . In the end, if no such 
inconsistencies are found, then our proposed explanation of the knee effect 
is ruled out. If evidence for missing energy beyond that expected in the 
standard model is found, then the cross section, its energy dependence and 
the energy dependence of the missing energy effect must be able to account 
for the knee. Of course it may be only part of the story behind the knee 
phenomenon, the rest coming from new astrophysics. 
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6 Summary and Conclusions 



We have employed several models, generic parameterizations, of new physics 
interactions to show that the knee in the cosmic ray spectrum in the region 
2 — 4 x 10 6 GeV can be described successfully with only a few parameters. 
We first wrote the total cross section as a sum of the contributions due 
to standard model and those due to new physics. In this simple picture, 
the observed showers are initiated either by the standard model interaction 
or by the new physics interaction, with weights given by the respective 
cross sections. We hypothesized that above the threshold for new physics a 
fraction y of the primary energy is in a form that the arrays do not detect, 
but that the observable energy is in showers similar to those modeled by the 
experiments. We fixed the spectral index of the incident flux at the value 
determined from the data below the knee. We found that the total proton- 
proton cross sections implied by the parameters determined by fitting the 
cosmic ray knee rise extremely rapidly in the energy range just below the 
knee. The result is that the discrepancy between the cross sections required 
to fit the knee with a single power law injection flux below about 10 9 GeV 
appear to be completely ruled out by the estimates of the cross section 
published by the Fly's Eye and AGASA collaborations in the 1980's and 
1990's. 

However the final pp cross section at ultra high energies is so large that 
the picture which leads to the simple form for the total cross section, Eq. [H 
is called into question. In this case multiple scattering is likely to dominate 
the nuclear collision and subsequent shower formation. In this limit, the 
nucleus becomes a black disc and one loses the clean separation between the 
new physics and standard model contributions and furthermore loses the 
direct link between the nucleon-nucleon and the nucleon-air cross section. 
Therefore it is not possible to test the new physics proposal purely on the 
basis of comparison with pp total cross section data. One would require a 
detailed analysis of the LHC data in order to determine if there exists a new 
interaction of sufficiently large strength along with sufficient missing energy 
in the final state to account for the knee in the CR spectrum. We pursued the 
consequences of this black disc picture for the Tibet III data. We modeled 
the new physics as invisible energy loss with a power law dependence on 
collision energy that sets in at a threshold energy Et, adopting a total cross 
section value extrapolated from our fit to accelerator data. Within this 
picture we also obtained a good fit to the cosmic ray flux in the vicinity of 
the knee. 

We conclude that, based on the published ultra- high energy total proton- 
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proton and proton-antiproton cross sections and a proper rescaling of the 
cosmic ray data, a consistent picture of the CR knee in terms of purely new 
physics effects with an invariant mass scale in the 1-3 TeV range emerges. 
Moreover, looking at the limit where high multiple scattering approaches a 
black disk picture, while the fraction of energy loss increases with energy, a 
natural solution to the knee in the CR spectrum again emerges. 
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